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Double quantum dots are convenient solid-state platforms to encode quantum information. Two- 
electron spin states can be conveniently detected and manipulated using strong quantum selection 
rules based on the Pauli exclusion principle, leading to the well-know Pauli spin blockade of electron 
transport for triplet states. Coherent spin states would be optimally preserved in an environment 
free of nuclear spins, which is achievable in silicon by isotopic purification. Here we report on a 
deliberately engineered, gate-defined silicon metal-oxide-semiconductor double quantum dot system. 
The electron occupancy of each dot and the inter-dot tunnel coupling are independently tunable by 
electrostatic gates. At weak inter-dot coupling we clearly observe Pauli spin blockade and measure a 
large intra-dot singlet-triplet splitting > 1 meV. The leakage current in spin blockade has a peculiar 
magnetic field dependence, unrelated to electron- nuclear effects and consistent with the effect of 
spin-flip cotunneling processes. The results obtained here provide excellent prospects for realizing 
singlet-triplet qubits in silicon. 



Gate-defined semiconductor quantum dots enable the 
confinement and manipulation of individual electrons and 
their spin [1 . Most of the relevant parameters - electron 
filling, energy splittings, spin states, exchange interaction 
- can be tuned in situ by electric and magnetic fields. Be- 
cause of this exquisite level of control, quantum dots are 
being investigated as candidate systems for spin-based 
quantum information processing [2] . In group III- V semi- 
conductors such as GaAs, the development of highly 
tunable double quantum dots has allowed the study of 
both single-electron and two-electron spin dynamics [3]- 
[?]• However, the nuclear spins always present in these 
materials produce strong decoherence of the electron spin 
degree of freedom and result in phase coherence times T2 
of below 1 ms [8^ 9^ . Conversely, group-IV semiconduc- 
tors such as silicon, silicon-germanium and carbon can be 
isotopically purified, leaving only spinless isotopes. The 
weak spin-orbit coupling [10 and the absence of piezo- 
electric electron-phonon coupling [11] allow for extremely 
long spin relaxation times Ti of order seconds, as already 
demonstrated in several experiments p!2HT4] . The phase 
coherence times have not been measured yet, but they 
are expected to reach ~ 1 s as well, in highly purified 
^^Si substrates with low background doping concentra- 
tion [15 . 

A widely successful method to observe and control spin 
phenomena in quantum dots [1 consists of defining a 
double quantum dot in a series configuration and tun- 
ing the potentials such that sequential electron transport 
requires a stage where two electrons must occupy the 
same dot. The eigenstates of a two-electron system are 
singlet and triplet spin states, separated by an energy 
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splitting AsT which can be large in tightly confined dots. 
The electron transport then becomes spin-dependent and 
can be blocked altogether when the two-electron system 
forms a triplet state [5,fT6l. This phenomenon, known 
as Pauli spin blockade, has been extensively exploited to 
investigate the coherence of single-spin [4 and two-spin 
states [3 in GaAs and InAs [17 quantum dots. There- 
fore, observing and controlling spin blockade in silicon 
is a key milestone to unravel the full potential of highly 
coherent spin qubits. Preliminary success has been ob- 
tained in Si [18 and SiGe [1^ devices, but in each case 
the double dot system under study resulted from local 
variations in the potential of a lithographically-defined 
single dot, making it difficult to control individual dot 
occupancies or inter-dot coupling. Spin-based quantum 
dot qubits require exquisite control of these parameters, 
so a highly tunable double-dot system in silicon is essen- 
tial. For singlet-triplet qubits in multivalley semiconduc- 
tors it is also crucial to ensure that a large valley-orbit 
splitting is present, to avoid the lifting of Pauli blockade 
due to valley degeneracy [20l |2TJ . 

Here we present an engineered silicon double quantum 
dot which shows excellent tunability and robust charge 
stability over a wide range of electron occupancy {m, 
n). The silicon metal-oxide-semiconductor (MOS) struc- 
ture utilizes an AI-AI2O3-AI multi-gate stack that enables 
very small dots to be defined, each with independent gate 
control, together with gate-tunable inter-dot coupling. 
Such multi-gate stacks have previously been used to con- 
struct single Si quantum dots with the ability to achieve 
single electron occupancy [22] . The double dot presented 
here exhibits spin blockade in the few-electron regime, 
from which we are able to extract a large singlet— triplet 
energy splitting and also investigate a new mechanism of 
singlet— triplet mixing in the weak-coupling regime. 
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FIG. 1. SEM and schematic view of the device, (a) 

Scanning electron micrograph of a device identical to that 
measured, (b) (Not to scale) Schematic cross-section view of 
the Si MOS double quantum dot. The architecture is defined 
by Bl, B2 and B3 (barrier gates), LI and L2 (lead gates), 
and PI and P2 (plunger gates). The gates are separated by 
an AI2O3 layer (light gray). Positive voltages applied to the 
lead and plunger gates induce an electron layer (black dashes) 
underneath the Si02. By tuning the barrier gates. Dot 1 and 
Dot 2 are formed. The coupling of the dots is adjusted using 
the middle barrier (B2). The regions coloured with red are 
the source (S) and drain (D) contacts formed via diffused 
phosphorus. 



I. RESULTS 

Device architecture. Figure [T] shows a scanning elec- 
tron micrograph (SEM) and cross-sectional schematic of 
the device, which incorporates 7 independently controlled 
aluminium gates. When a positive bias is applied to the 
lead gates (LI and L2) an accumulation layer of electrons 
is induced under the thin Si02, to form the source and 
drain reservoirs for the double dot system. A positive 
voltage on the plunger gate PI (P2) causes electrons to 
accumulate in Dot 1 (Dot 2). Independent biasing of PI 
and P2 provides direct control of the double-dot electron 
occupancy {m, n). The tunnel barriers between the two 
dots and the reservoirs are controlled using the barrier 
gates: Bl, B2 and B3. The middle barrier gate B2 deter- 
mines the inter-dot tunnel coupling. The electrochemical 
potentials of the coupled dots can also be easily tuned to 
be in resonance with those of the source and drain reser- 
voirs. As shown in Fig. [ijb), gates LI and L2 extend 
over the source and drain n+ contacts, and also overlap 
gates Bl and B3. The upper-layer gates (PI and P2) 
are patterned on top of the lead and barrier gates. The 
lithographic size of the dots is defined by the distance 
between adjacent barrier gates (~30 nm) and the width 
of the plunger gates (^50 nm), as shown in Fig. [TJa). 

Inter-dot tunnel coupling tunability. Figure |2] 
shows the measured differential conductance of the de- 
vice as a function of the plunger gate voltages, Vpi and 
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FIG. 2. Characteristics at different inter-dot tunnel 
coupling. Measured stability diagrams and energy landscape 
of the double dot system ranging from weak to strong inter- 
dot tunnel coupling (a) — (c) and (d) — (f) respectively, for Vli 
= Vi2 = 3.0 V, V^Bi = 0.76 V, Vbs = 1.0 V and Fsd = 0. 
From lower to higher Vb2, the tunnel barrier height decreases 
resulting in stronger inter-dot tunnel coupling, (a) A checker 
box pattern, (b) honeycomb pattern and (c) diagonal parallel 
lines indicate that the two dots merge into a single dot as the 
coupling is increased |23j. 



Fp2, with all other gate voltages held constant, together 
with sketches of the energy landscape of the double dot. 
The charge-stability maps moving from Fig.|2];a) to[2|c) 
clearly show the effects of an increasing inter-dot coupling 
as the middle barrier-gate voltage Fb2 is increased, lower- 
ing the tunnel barrier between the dots. Fig. [2|b) shows 
the characteristic honeycomb-shaped stability map rep- 
resenting intermediate inter-dot coupling [23], obtained 
at Fb2 = 1-32 V. At lower middle barrier-gate voltage, 
Fb2 = 1-20 V, we observe a checker-box shaped map 
[Fig. |2|a)], since the middle barrier is opaque enough 
to almost completely decouple the two dots. In contrast, 
the stability map in Fig. [2|c) shows the formation of di- 
agonal parallel lines at Fb2 = 1-40 V. Here the two dots 
effectively merge into a single dot due to the lowering of 
the middle barrier [Fig. [2|f)]. The transport measure- 
ments shown here do not allow a precise determination 
of the electron occupancy {m, n) in the dots, since it 
is possible that electrons remain in the dots even when 
/sD is immeasurably small. For the regime plotted in 
Fig. [2] there were at least 10 electrons in each dot, based 
on our measurement of Coulomb peaks as we further 
depleted the system. An absolute measurement of dot 
occupancy would require integration of a charge sensor 
into the system [7]. These results nevertheless demon- 
strate that the multi-gated structure provides excellent 
tunability of coupling while maintaining charge stability 
over a wide range of electron occupancy. 

Capacitances and charging energies. Application 
of a DC source-drain bias Vsd causes the triple-points 
in the weakly-coupled regime [Fig. |2|a)] to extend to 
form triangular shaped conducting regions [Fig. [sj^a)] 



3 



On 



(a) 

2.124- 
2.104- 



/sd(pA) OE 



/sd(pA) 



5 = 0T 
v' FsD=-2.5mV 



2.140 
> 

2.120 



= + 2.5 mV 
B = OT 




2.008 2.024 

^P2 (V) 



FIG. 3. Pauli spin blockade at weakly coupled regime. 

Current /sd as a function of Vpi and Vp2 for B = T. The 
lead and barrier gate voltages were fixed at Vli = Vl2 = 3.2 V, 
Vbi = 0.656 V, Vb2 = 1.176 V and Vb3 = 0.940 V throughout 
the experiment, (a) For Vsd = — 2.5 mV, the ground state and 
excited states of a full bias triangle are shown. The current 
flows freely at the S(0,2)— S(l,l) transition as illustrated in 
the box marked by red dot. (b) The same configuration at 
Vsd = +2.5 mV, the current between the singlet and triplet 
states is fully suppressed by spin blockade (green star box) 
except on the bottom (blue cross box) of the bias triangle. 
The blue cross box shows how a leakage current arises the 
Pauli spin blockade region. 
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FIG. 4. Singlet-triplet splitting, (a) — (c) DC measure- 
ments of the triangle pair analysed in Fig. [3] at Vsd = 
+2.5 mV, for different in-plane magnetic fields, B (scale bar 
same as Fig. |3jb)). The singlet— triplet splitting, Ast, is de- 
fined by the triplet and singlet state of (0,2) as depicted in 
(a). As the magnetic field increases, Ast decreases along the 
detuning axis of the triangle [labeled e in (b)]. (d) The energy 
spacing Ast as a function of in-plane magnetic field B. Ast 
decreases at a rate ^0.12 meV/T and is expected to approach 
zero at 11.3 T. From the linear fit (red line) through Ast, the 
^-factor is 2.1 ± 0.2. 



from which the energy scales of the double dot sys- 
tem can be determined |23J. From a triangle pair, we 
extract the conversion factors between the gate volt- 
ages and energy to be ai = eFso/^^'^Pi = 0.089e and 
a2 = eVsB/^Vp2 = 0.132e, where SVpi and SVp2 are the 
lengths of the triangle edges, as shown in Fig. |3|a). The 
total capacitances of Dot-1 and Dot-2 can then be calcu- 
lated [23 , giving Ci = 16.3 aF and C2 = 14.5 aF. The 
charging energies of the two dots are then Eq^i = e^/Ci 
= 9.8 meV and £^c,2 = e^/C2 = 11 meV, indicating that 
the left dot is slightly larger than the right dot. 

Pauli spin blockade. Figure [3] shows the current 
/sD through the double dot as a function of the two 
plunger gate voltages when measured with both posi- 
tive [Fig. [3]^ a)] and negative [Fig. |3jb)] source-drain bi- 
ases. Here we observe a suppression of current at one 
bias polarity, the characteristic signature of Pauli spin 
blockade [H [7]. At Vsd = —2.5 mV we observe a pair 
of overlapping full bias triangles, as shown in Fig. |3ja). 
Resonant transport through the ground state and the 
excited states in the double dot occurs when the states 
within the dots are exactly aligned, leading to peaks in 
the current which appear as straight lines parallel to the 
triangle base in Fig.jsja). The non-resonant background 
current level at the centre of the triangle is attributed 



to inelastic tunneling. The non-zero current throughout 
the triangular region indicates that electrons from the 
reservoir can tunnel freely from the S(0,2) singlet state 
to the S(l,l) singlet state, as depicted in the cartoon (red 
box in Fig. [3|. Note that here we define {m, n) as the 
effective electron occupancy [18 , while the true electron 
occupancy is (m+mo, n+no). The Pauli blockade ex- 
pected for two-electron singlet and triplet states occurs 
when the total electron spin of each dot is zero in the 
(mo, no) state. 

At the complementary positive bias of Vsb — +2.5 
mV we observe strong current suppression in the region 
bounded by the dashed lines in Fig. [sj^b). The suppres- 
sion arises because the transition from T(l,l) to S(0,2) is 
forbidden by spin conservation during electron tunneling. 
Once the T(l,l) triplet state is occupied, further current 
flow is blocked until the electron spin on the right dot re- 
verses its orientation via a relaxation process (green star 
box in Fig.[3| [6 , 7 . 

Singlet-triplet splitting. In a magnetic field B there 
are four accessible spin states: the singlet S; and three 
triplets T_, To and T+, corresponding to 5'z = —1, 0, 
+1. The singlet— triplet splitting Ast is the energy dif- 
ference between the blockaded ground state S(0,2) and 
the excited state T_(0,2) [7l[T8]. Here we study Ast as 
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a function of 5, applied parallel to the substrate, by mea- 
suring spin blockade at a positive bias. Figures [4]^a— c) 
show the bias triangles in the spin blockade regime at 
increasing magnetic fields 5 = 2, 4 and 6 T, with the 
splitting AsT marked in Fig. [4]^a). The measured split- 
ting AsT decreases linearly with increasing B [Fig.[4]^d)], 
as expected, since the triplet states split linearly by the 
Zeeman energy, Ez = ±5'z|^|/iB^, where /i^ is the Bohr 
magneton and iSz is — 1, 0, +1. A linear fit through 
^st{B) yields a Lande ^-factor of 2.1 ± 0.2, consistent 
with electrons in silicon. 

We observe an exceptionally large value of the (0,2) 
singlet-triplet splitting at 5 = 0, Ast ^ 1-4 meV. 
This result is striking because it implies that the near- 
est valley-orbit state must be at least 1.4 meV above the 
ground state. The first excited valley-orbit state should 
be a combination of the ±z valleys. It would lift the spin 
blockade [20, ^l^^ and show no remarkable energy shift 
in a magnetic field, in contrast with our observations. 
Therefore, such a state must lie above the triplet state 
we observe in Fig.[4j The ability of our structures to gen- 
erate such a large valley-orbit splitting removes a major 
concern on the realizability of singlet-triplet qubits in a 
multivalley material such as silicon. 

Leakage current in blockade regime. If some 
mechanism exists to mix the singlet and triplet states 
or to induce transitions between them, then the spin 
blockade can be lifted, leading to a measurable leakage 
current [6^. As shown in the blue cross box in Figure [sj 
transitions from T(l,l) to S(l,l) can lift the blockade, al- 
lowing electrons to transit the double dot until the next 
triplet is loaded, resulting in a non-zero time-averaged 
leakage current. Fig. [sja) shows the surface plot of the 
leakage current /sd as a function of both detuning e and 
magnetic field 5, while Figs, jsjb) and [sjc) show line 
traces of /sd as a function of B at zero detuning and /sd 
as a function of e at zero magnetic field, respectively. We 
find that the leakage current has a maximum at 5 ~ 
and falls to zero at |5| - 700 mT. 



II. DISCUSSION 

The suppression of leakage current by an applied mag- 
netic field has been observed in GaAs double quantum 
dots (DQDs) [6 and attributed to the effect of hyperfine 
coupling between the electron spins and the surround- 
ing bath of nuclear spins. In that case the width SB of 
^sd(^) yields the average strength of the hyperfine field. 
For an unpolarized nuclear spin bath 5B ^ (^^max/V^, 
where ^^max is the hyperfine field assuming fully polar- 
ized nuclei and N is the number of nuclei overlapping 
with the electron wave function. For a typical GaAs 
dot overlapping with ^ 10^ — 10^ nuclei, SBmax ~ 6 T 
=^ 6B 2 — 6 mT ^ [71 [24]. In natural silicon, however, 
the hyperfine interaction is much smaller than in GaAs, 
with ^^max ^1.9 mT [25]. Therefore, hyperfine coupling 
can be ruled out as a mechanism for the leakage current 
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FIG. 5. Spin-flip cotunneling in Pauli spin blockade 
regime, (a) A surface plot of leakage current through spin 
blockade as a function of energy detuning e and magnetic field 
with gates settings as in Fig[3|;b). (b) Cut along B at e = 
energy detuning axis (black arrow) while (c) Cut along e at 
B — field (blue arrow) . Fits of experimental data with spin- 
flip cotunneling relation give Fd = 34 /xeV, t = 0.5 /xeV and 
T 115 mK. 



here. 

An alternative mechanism for a transition from triplet 
to singlet has been recently proposed, where the spin flip 
is caused by inelastic cotunneling [26] . The spin-fiip rates 
due to cotunneling from the spin-polarized triplet states, 
T±(l, 1), are exponentially suppressed when the Zeeman 
energy is large compared to the thermal broadening of 
the electron states in the leads (i.e., for gji^B > k^T^ 
where T is the electron temperature and B is the applied 
magnetic field). A rate-equation analysis accounting for 
the energy dependence of the spin-fiip cotunneling rates 
[26] then gives a simple form in the limit of weak inter- 
dot tunneling t and weak cotunneling W^^^ compared to 
the tunnel rates rs,D between a dot and its nearby source 
or drain lead {V^t < k^T, W^^^ <C rs,D) [22]: 



/sD(5) = e, 



3Wc'ot/2 



sinh 



( gfiBB \ 



, e = 0. (1) 



Here, the B = spin-fiip cotunneling rate (for kBT > 
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V2i and |e| < |A|, |e|FsD) is: 
kBT 



A 



A - 217' - 2\e\VsD 



(2) 



with mutual (inter-dot) charging energy and A = 
aiSVpi + a2SVp2 for plunger gate voltages (^Ii5i,p2 mea- 
sured from the effective (0,1) — (1,1) — (0,2) triple 
point (lower-left corner of the bias triangle in Fig.jsjb)). 
Eq. ([2| accounts for virtual transitions between effective 

(1.1) and (0, 1) (first term) as well as effective (1,1) and 

(1.2) charge states (second term). 

In the present case, A |e|VsD ^ U\ The higher 
current level in the upper right corner of Fig. [sj^b) fur- 
ther suggests Fd ^ Fs, giving (for this particular exper- 
iment): 



cot 



TrTi \|e|VsD 



(3) 



Using the above expression for then use Eq. 

to fit to the Isb{B) data in Fig. pjb), giving us Fd = 
34 /ieV for the tunneling rate and T = 115 mK for the 
electron temperature. 

The B = spin-flip cotunneling rate W^^^ is energy- 
independent in the limit < k^T. However, the leak- 
age current does acquire a dependence on the energy de- 
tuning, e = aiFpi — OL2Vv2^ when the escape rate from 
the double-dot due to resonant tunneling is suppressed 
below the spin-flip cotunneling rate. This leads to a 
Forentzian dependence of the current on detuning e with 
a t-dependent width ^e: 



/sD(e) 



^c'ot 



31 

SFd^ 

^c'ot 



(e/Je)^ 

.2\ 1/2 



5 = 0, 



(4) 
(5) 



Eq. Q is valid in the same limit {^t < k^T, W^^^ <C 
rs,D) as Eq. (Q). In the strong-tunneling limit, 
V^t > /cbT, the theory predicts that /(e) should show 
a strong resonant-tunneling peak of width ~ t, followed 
by a slowly- varying Forentzian background described by 
Eq. (Q) at large e. The absence of a strong resonant- 
tunneling peak in the data of Fig. Isjc) confirms that the 
device is operating in the regime ylt < k^T^ justifying 
our use of Eqs. ([T]) and Q to analyse the data. 

A nonlinear fit to the /sD(e) data in Figs.jsjc) using Eq. 
Q yields t = 0.5 /ieV for the inter-dot tunneling rate, 
using our previously determined values Fd = 34 /ieV and 
T = 115 mK. These parameter values are well within the 
experimentally expected range. The small value of t in- 
dicates weak inter-dot tunnel coupling, consistent with 
the results shown in Fig. [sj^b). We conclude that the 
spin-flip cotunneling mechanism provides a fully consis- 
tent explanation of the observed leakage current in the 
spin blockade regime. The mechanism could be applied 



to reanalyse previous experiments in group IV semicon- 
ductors [28j where the nature of the leakage current was 
not fully understood. 

In conclusion, we have presented a lithographically- 
defined double quantum dot in intrinsic silicon showing 
excellent charge stability and low disorder. The multi- 
gate architecture provides independent control of elec- 
tron number in each dot as well as a tunable tunnel 
coupling. We observed Pauli spin blockade in an effec- 
tive two-electron system from which we extracted the 
singlet— triplet splitting. The leakage current in the spin 
blockade regime is well explained by a spin-flip cotunnel- 
ing mechanism, which could be of widespread importance 
in group-IV materials with weak hyperfine coupling. The 
results obtained here provide a pathway towards inves- 
tigation of spin blockade in silicon double quantum dots 
with true (1,1) and (2,0) electron states. Towards this 
end, we are planning future experiments incorporating 
a charge sensor to monitor the last few electrons j29] . 
We anticipate that such an architecture will provide ex- 
cellent prospects for realizing singlet— triplet qubits in 
silicon [30j . 



III. METHODS 



Fabrication steps. The devices investigated in this 
work were fabricated on a 10 kl]-cm n— type high resistiv- 
ity (100) silicon wafer using standard micro- fabrication 
techniques. The n+ source and drain ohmic contacts 
regions in Fig. [TJb) were produced via high concen- 
tration phosphorus diffusion at ~ 1000 °C, resulting in 
peak dopant densities of ~ 10^^ cm~^. Next, the high- 
quality Si02 of 10 nm thickness was grown via dry ther- 
mal oxidation in the central region at 800 °C in O2 and 
dichloroethylene. The barrier gates were first patterned 
on the thin Si02 region using electron beam lithography 
(EBF) followed by thermal evaporation of 40 nm thick 
aluminium and lift-off process. Before the next EBF step, 
the barrier gates are exposed to air for 10 mins at 150 °C 
to form ~4 nm of AI2O3 acting as a dielectric layer. This 
process was repeated for lead gates and plunger gates 
layers with aluminium thicknesses of 40 nm and 120 nm 
respectively. A final forming gas anneal (95% N2 and 
5% H2) was performed for 15 mins to achieve a low den- 
sity of Si-Si02 interface traps, of order 10^^ cm~^ eV~^, 
as measured on a similarly processed chip [31]. The low 
trap density is clearly reflected in the device stability and 
the low level of disorder observed in the transport data 
shown in the results section. 

Experimental setup. Electrical transport measure- 
ments were carried out in a dilution refrigerator with a 
base temperature T ~ 40 mK. We simultaneously mea- 
sured both the DC current and the differential conduc- 
tance dl/dVsd^ the latter using a source-drain AC exci- 
tation voltage of 100 jaV at 87 Hz. 
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